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Abstract. We consider a class of fourth order elliptic systems which include the Euler- 
Lagrange equations of biharmonic mappings in dimension 4 and we prove that weak limit 
of weak solutions to such systems is again a weak solution to a limit system. 



1. Introduction 

In this paper we consider fourth order elliptic systems of equations of the form 

(1.1) A 2 u = A(D ■ Vu) + div(E ■ Vu) + G ■ Vu + AVt ■ Vu 
for an unknown map u: B 4 C M 4 — > W 71 , i.e., in components, 

(1.2) A V = A (pid a J) + d a (E%dpu>) + G%d a u* + A{tf°)d a ui 

where a, /3 — 1, 2, 3, 4 and i,j — l,...,m. The coefficient functions D, E, G, Q are assumed 
to satisfy 

(1.3) DeW 1 ' 2 , E G L 2 , GGL 4/3 , fieW /1 ' 2 (B 4 ;so(m)®A 1 M 4 ). 

We study compactness of the space of solutions in the weak sequential topology of the 
Sobolev space W 2,2 . 

Let us note immediately that under the above assumptions G ■ Vw is just in L 1 , as 
Vu G W 1,2 C L 4 by Sobolev imbedding in dimension 4. Thus, (II. ip is critical. However, 
the whole point is that we allow D, E, 67, Q to depend nonlinearly on u. The class of 
systems we consider contains, in particular, the Euler-Lagrange equations of biharmonic 
maps from domains in M 4 into compact Riemannian manifolds. Our approach relies in a 
crucial way on the antisymmetry of the 1-form Q and on the use of nonlinear counterparts of 
the Hodge decomposition, originating in gauge theory. This key idea is due to T. Riviere; it 
has been first used in his pioneering paper [7] on conformally invariant second order elliptic 
systems in the plane, with harmonic maps from planar domains into compact Riemannian 
manifolds serving as the crucial example. Later on, Riviere-Struwe [8], Lamm- Riviere [3], 
and Struwe [10] extended this approach to stationary harmonic maps in higher dimensions 
and to biharmonic maps. 



Date: January 8, 2013. 

All authors partially supported by MNiSW Grant no N N201 397737, Nonlinear partial differential 
equations: geometric and variational problems. 

1 



2 PAWEL GOLDSTEIN, PAWEL STRZELECKI, AND ANNA ZATORSKA-GOLDSTEIN 

Let us also note that ( II. ip in full generality is wider than the class of biharmonic maps. 
It can happen in dimension 4 that a solution of ( II .ip is continuous, even C 1 ^, but still not 
C 2 . We shall comment on that later on; let us now state the main result. 

Theorem 1.1. Suppose (u k ) is a sequence in W 2,2 (M 4: , IR m ) of weak solutions to the system 

(1.4) A 2 u k = A(D k ■ Vu k ) + div(E k ■ Vu k ) + G k ■ Vu k + AQ k ■ Vu k in B 4 . 

Suppose u k — ^ u in W 2,2 (M 4 , M m ). If the coefficients D k , E k; G k , Q k are weakly convergent 
in their respective Sobolev spaces, i.e. 

(1.5) D k ^D, Q k ^Q, E k E, G k ^ G, 
then u is a weak solution to the limit system 

(1.6) A 2 u = A(D ■ Vu) + div(E ■ Vu) + G ■ Vu + AQ ■ Vu in B 4 . 

Related compactness results for biharmonic maps, along with an energy identity account- 
ing for the possible 'energy loss' under the limit passage have been obtained by Hornung and 
Moser [2] and Laurain and Riviere jl]. (For second order elliptic systems Au k = Q k ■ Vu k , 
Sharp and Topping [9] provide a compactness theorem under an assumption that rules out 
the concentration of the energy J |Vufc| 2 but allows for concentration of J |V 2 w*;|.) 

One of our main points is that the proof in this paper is different from [2] and jl]: 

contrary to these two papers, we do not rely at all either on continuity of solutions (or 

their first and second order derivatives) or on any other improved regularity properties (like 

higher integrability of V 2 u), working all the time just in W 2,2 and using the imbedding 
W 2,2 c jyM.rj 

Besides K. Uhlenbeck's Theorem I2.1| the main tool is the concentration- 
compactness method of P.-L. Lions. The combination of the two allows us, very roughly 
speaking, to reduce the complexity of (jl.4p ~ (jl.6p to the case 

A 2 u k = G k ■ Vu k , u k -± u in W 2 ' 2 , 

where G k is a bounded sequence in L 4//3 , with G k — ^ G. Then, a passage to the limit (in 
dimension 4) can be justified by an application of Sobolev inequality, the concentration- 
compactness lemma j5j Lemma 1.2], and a standard capacity type argument. 

Remark 1.2. Conditions (II. 3p imposed on the coefficients allow for nonlinear dependence 
on u. Roughly speaking, the class ( II. ip contains systems of the form A 2 u = R(u, Vw, V 2 w) 
where the right hand side R depends polynomially on (the entries of) Vii and V 2 m, with 
coefficients that are smooth and bounded in u, so that 

\R(u, Vw,V 2 w)| < |V 2 w| 2 + |Vw| 4 + l. 

The point is that the terms of R depending quadratically on V 2 n need to have some 
structure, whereas all the other terms are allowed to be arbitrary. A model case looks as 



In Section 2] we provide an example showing that admits weak solutions that are C 1 but not C 2 . 
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follows (summation over repeated indices is understood): 

(1.7) D = D(u, Vu) G W 1 ' 2 , D% = d?(u)d tt u l ; 

(1.8) E = E{u, Vu, V 2 u) G L 2 , E% = ef'^ {u)d a {d p u l ) + e%®{u)d a u%u a ; 

(1.9) G = G(u, Vu, V 2 u) G L 4/3 , Gl = gl p>M {u) (d^u 1 ) + dpu l d,u s )d a u p ; 
finally, the antisymmetric term Q = Q(u, Vu) G H /1,2 (B 4 ; so(m) Cg> A 1 !! 4 ) is given by 

(1.10) fig = Wl d a w j - wfdnW* = -n£ . 

where w = v o u of u for some bounded smooth map v. M. m — > R m . 

All the coefficients d,e,g (with various indices) in (jl.7D — (jl.9j) are assumed to be of class 
C 1 fl W 1,oc on M. m . Struwe [T0| Section 2] explains that the biharmonic map equation can 
be written in that form, with v being the normal to the target manifold. In that case, 
d, e, g depend explicitly on v and its derivatives, and the growth estimates 

|D| + |fi|<|V«|, 

(1.11) \E\ + \VD\ + |VO| < |V 2 n| + |Vn| 2 , 

\G\ < |V 2 n||Vn| + |Vn| 3 

follow from (jl.7p — (ll.lOp . Under these assumptions, for every weakly convergent sequence 
Uk — ^ u in W 2 ' 2 we have convergence (jl.5p . with 

L 2 3 E k = E{u k , Vu k , V 2 u k ) ^E = E{u, Vu, V 2 u) 

in L p for p < 2, and L 4 / 3 3 G k = G(u k , Vu k , V 2 u k ) -± G = G(u, Vu, V 2 u) in L s for 
s < 4/3. Since E k is bounded in L 2 and G k is bounded in L 4 / 3 , and we deal with a bounded 
domain, it is an exercise to see that in fact E k — E and G k — ^ G also for the limiting 
exponents p = 2 and s = 4/3. 

We do not need the full strength of (O HOI) . 

Remark 1.3. Our proof depends in a crucial way on Sobolev imbedding in dimension 4. 
It would be interesting to know what happens in higher dimensions. For example: is the 
convergence u k — > u in BMO, combined with the boundedness of (u k ) in W^ 2,2 , sufficient 
to guarantee (a) that u solves the limiting system, (b) that u k — >• u strongly in W 2,2 7 It 
is possible to check, using the sharp version of Gagliardo-Nirenberg inequality ||Vio|| 2 4 < 
H'HIbmoII'MI vk 2 > 2 (°f- Meyer-Riviere |6j, or [11]). that both answers would be positive for 
uniformly bounded weak solutions of the simplified system A 2 u k = G k Vu k , with G k = 
G(u k ,Vu k ,V 2 u k ), where G is of the form (II. 9p . It seems plausible that convergence of, 
say, biharmonic maps in BMO prevents bubbling and loss of energy in the limit. 

The rest of the paper is organized as follows. In Section 2, we recall Uhlenbeck's result 
in the form that is used later on. Section 3 forms the bulk of the paper. There, we first 
explain the strategy of the proof in more detail, and then carry out the necessary estimates, 
pass to the limit and remove the singularities of the limiting system. Finally, in Section 4 
we give an example showing that (11. ip can have solutions in C 1 \C 2 . 
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2. Uhlenbeck's result 

We remind that in our case we consider a mapping u going from a ball B 4 C M 4 into 
M. m . Below we state the theorem of Uhlenbeck in a form adjusted to the situation. Note 
that if ft G W 1>2 (B 4 ) then ft G L 4 (B 4 ) because of the Sobolev imbedding theorem. 

Theorem 2.1. There exist a number e = e(m) > and a constant C > such that for 
any ft G W 1,2 (! 4 ; so{m) <g> A X R 4 ) which satisfies 

||Vft|| L 2 + ||ft||l4 < e 
there exist P G W 2 ' 2 (M A , SO(m)) and £ G I¥ 2 ' 2 (B 4 , so{m) <g> A 2 M 4 ) such that 

(2.1) (dP) P" 1 + PUP" 1 = *di on B 4 

and d(*f) = on B 4 , £| aB 4 = 0. 
Moreover 

(2.2) ||V 2 P|| L2 + ||VP|| X 4 + ||V 2 £|b + ||V^|| L 4 < ||Vft|| L2 + ||ft|| 2 4. 

Uhlenbeck's Theorem is, in fact, a local theorem in the sense that we can use it not 
only on the unit ball B 4 , but on any ball, and, as long as we consider balls with uniformly 
bounded radii, we can choose the constant e in an uniform way (i.e. independently on the 
radius of the ball). This is in accordance with the original use of this theorem to prove the 
existence of global Coulomb gauges on compact manifolds. Indeed, a look at the proof of 
Lemma 2.5 in [13] shows that we can choose e = (2(Cpr + l)Cs')~ 1 , where Cp and Cs are 
the constant in the Poincare and Sobolev inequalities for the unit ball, and r denotes the 
radius of the ball B. Thus e = (2(CpR + l)^) -1 can be chosen as a uniform estimate for 
all balls with radius bounded by R. 

Corollary 2.2. Theorem \2.1\ holds for any ball B C B 4 in place o/B 4 , and the constant e 
can be chosen uniformly for all such balls. 

Another corollary deals with the problem of weak continuity of P and £ with respect to 
ft. Note that Theorem 12.11 does not claim that either of them is defined uniquely. 

Corollary 2.3. Suppose {ft^} is a sequence in W 1,2 (M 4: ; so(m) <8> A 1 M 4 ). Assume 

(2.3) ft fc ft in W 1 ' 2 . 

Assume that Pk and are chosen so that (12. ip and (12. 2p of Theorem \2.1\ hold with ft^ in 
place o/ft. 

Then both Pk and ^k are uniformly bounded in W 2 ' 2 , and for any subsequence o/(ftfc) for 
which Pk and are weakly convergent in W 2,2 to P and £, respectively, conditions (12. ip 
and (12. 2p of Uhlenbeck's Theorem do hold for ft, P and £. 

In other words, the decomposition for the limit matrix ft can be effectuated with (any) 
weak limit of the transformations Pk and forms 

Indeed, assume ft^ — ft w W 1,2 and that ft& satisfy, uniformly with respect to k, 
assumptions of Theorem 12.11 We obtain then, by the theorem, transformations Pk and 
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such that 

(2.4) dP k + P k n k = *d& P k on B 4 
and 

(2.5) ||v 2 p fc || L2 + \\vp k \\ L 4 + \\v%\\ L 2 + ||v&||£4 < \\vn k \\ L2 + \\n k \\ 2 L4 < m, 

where M > is a constant which is independent of k. 

Since, by (I2.5p . the sequences (£ k ) and P k are bounded in W 2 ' 2 , we can choose subse- 
quences (for simplicity still indexed by k) such that 

£fc £ weakly in W 2,2 , 

P k — ^ P weakly in W 2,2 and a.e. 

Thus (after again choosing subsequence) we may assume that 

d£ k — > d£ strongly in L 4_<5 , 

dP k — > dP strongly in L 4_<5 , 

£l k — > £1 strongly in L A ~~ 5 . 

for any small 5 > 0. Since P k are also uniformly bounded in L°°, we can take the L 4_<5 -limit 
on both sides of (12. 4p obtaining 

(2.6) dP + Ptt = *di on B 4 

in the sense of distributions, which proves that P, £ and Q indeed satisfy (12. ip . The 
remaining estimates and boundary conditions in Theorem 12.11 are obvious. 

This sort of continuity of the decomposition of Q, i.e. the fact that the transformation 
of the limit Q may be attained by taking the weak limits of the elements of decomposition 
of Q k allows us later to estimate the W 2,2 -norm of differences between the elements of 
decomposition of Q k and Q. 

3. Proof 

Let us first give a rough sketch and plan of the proof. The key idea is to prove that u 
solves the limiting system (11.61) outside a countable set of points and then remove these 
possible singularities with the use of a properly chosen test function. A standard argument, 
cf. Section 13. 1\ shows that B 4 \ Ai, where Ai is finite, can be covered by balls Bj such 
that ll^jfellw 1 . 2 ^-) is small, so that Uhlenbeck's decomposition can be applied inside each 
Bj separately. Next, in Section 13.21 we fix one of the Bj and, following the crucial ideas 
of Lamm-Riviere [3] and Struwe [10| . we use the equation (dP k )P k ~ l + PkVLkP^ 1 = *d£ k on 
Bj to rewrite (jl.4p as 

(3.1) A(P k Au k ) + a perturbation = H k ■ Vu k + *dA£ k ■ (P k Vu k ) on Bj, 

where H k depends on D k , E k ,G k , P k , Q k and their derivatives. 

It is a purely routine matter to pass to the weak limit on the left hand side of (13. ip . On 
the right hand side, after passing to subsequences, H k is bounded in L 4 / 3 and converges 
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strongly in all L s for s < |, and is bounded in L 4 by the Sobolev imbedding, and 
converges strongly in L q for all q < 4. Thus, H k • Vitfe is bounded in L 1 but a priori does 
not have to converge in T>' . The second term, *dA£/t ■ (PfcVitfc), presents a similar difficulty 
To cope with that, we apply in Section |3~31 P.-L. Lions' concentrated compactness method 
[5], following earlier work of Freire, Muller and Struwe [T| on wave maps and harmonic 
maps, and the second and third author [12] on if-systems. The key idea is to exploit the 
existence of second order derivatives of u. This yields 

H k ■ Vu k + *dA£ k ■ (P k Vu k ) — >H -Vu + *dA£ ■ (PVw) + Sj on Bj, 

where H ■ Vu + *dA£ ■ (PVw) is the desired term of the limit system (11.6P rewritten in the 
(P, £) gauge, and Sj is a combination of Dirac delta measures, supported on a countable 
subset A of By 

To complete the proof, in Section I3~41 we show that each Sj must be zero, using a capacity 
argument, based on the fact that W 2 ' 2 (M. A ) contains unbounded functions. Thus, the limit 
u of Uk satisfies (I1.6P in B 4 \ A\\ another application of the same argument shows that A\ 
must be empty. 

3.1. Preparation to Uhlenbeck's transformation. Since the sequence is weakly 
convergent, it is bounded in W 1,2 ; we shall denote the bound on its norm by M. 

In what follows, we want to cover the ball B 4 by balls Bj in such a way, that in every of 
the Bj we may, after passing to a subsequence, assume that ||fifc||w , i,2(£ j ) < e, where e is as 
in Corollary !2.2l to Uhlenbeck's Theorem. This might not be possible for the whole B 4 , but 
it is outside a finite set of points. To visualize this better, replace B 4 by a four-dimensional 
cube and consider its dyadic decomposition into cubes Cij, where the second subscript j 
counts subsequent cubes of a specified generation i. 

A cube C is bad, if one cannot choose a subsequence of such that its H /1,2 -norm is 
bounded on C by e, i.e. if, for all k sufficiently large, ||fifc||w'i.2(c) > £ - Notice that, in 
every generation of the dyadic decomposition, the number of bad cubes C^j is bounded by 
the same constant iV = [Me -1 ], and that if a cube is not bad (i.e. it is good), neither are 
all its descendants. The intersection of all bad cubes, 

(3.2) A 1 :=f}( \J bad cubes C it j\ 

i=l ^l<j<N ' 

is a finite set of (at most N) points, and any point not in A\ lies in a cube that is good. 

These points are, in fact, accumulation points of the W^^-norm of the weakly convergent 
sequence (Qk), more precisely — these accumulation points, for which the energy loss 
exceeds e. 

For the next two subsections, let us fix an arbitrary good ball B = Bj contained in 
B 4 \ A\. We consider (not changing the notation, to keep things simple), instead of the 
whole sequence u^, only the subsequence of (and of D^, E^, and Q^) f° r which 
H^fcllw^CB-) — £■ By Corollary 12.21 we can assume that the Uhlenbeck's theorem holds for 
B, so there exist and such that dPk P k _1 + PfcfifcP^T 1 = *d£k- 



WEAK COMPACTNESS FOR 4-TH ORDER SYSTEMS 7 

3.2. Transformation of the equation: calculations. The calculations below follow 
closely and with more detail the brief calculations by Struwe in [TU]. We provide them 
for the reader's convenience, and also because we shall need some more knowledge on the 
structure of certain terms in our reasoning. 

We recall the indices of the multidimensional objects that appear in our system: 



a 1 1 



u = (u i ), Vu = {d a u i ), P={PV), D = {D, 

e = (e%), G = Gt n = (n«), 

with i,j = l,...,n and a, (3 = 1,2,3,4. To simplify the notation without making the 
calculations ambiguous we shall use the standard summation conventions. 

Furthermore, it is often convenient to omit at least some of the indices. In that case, 

• multiplication of tensor objects that is denoted by a dot (•) is a standard scalar 
product in the M. 4 x W m space, e.g. 



E ■ Vu = (E%d,u% = £ E % d ?A 

V /3=l 3=1 J 

(An v i 



""=1 3=1 

multiplication of tensor objects that is not denoted by any operator sign is standard 
matrix multiplication, e.g. 



PA 2 u = (P ij A 2 u j ). = (j2P ij A 2 uA 



tensor multiplication ((g)) denotes tensor product in M 4 (and then, possibly, matrix 
multiplication in the coordinates), e.g. 



VP (g Vu = VP ij (g Vu j = ( y (d a P ij ) a (g (d p v?) i 



j=l / 

Below, we transform the system (jl.4p for Uk- In the calculations that follow we omit the 
index k (one should not confuse this temporary notational simplification with the claim 
that u - a weak limit of (uk) - satisfies (jl.6P ; proving this is the goal of our paper). 

Applying P to the Laplacian on the left-hand side of the system, we obtain 

A(PAuY = A{P ij Au j ) = d a d a {P l] dpd p u j ) 

= p ij d a d a dpdpu j + dpd a d a p ij dpu j - dp(d a d a p ij dpu j ) 

+ 2d a d (d a P ij dpu j ) - 2d a (df,d tt F i 'dffu>), 
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which can be rewritten shortly as 

A(PAu) 1 = (PA\) ! + (V(AP) ■ V«)' - div(APVu) 1 

+ 2div 2 (VP <g> Vu) 1 - 2div(V 2 P ■ Vu)\ 

By the equation (11.41) satisfied by Uk we have, still omitting the index k for the sake of 
simplicity, 

(PA 2 «) ! = P i] A(D ■ Vu) j + P i Miv(P • Vu) j + P' ll {G ■ Vu) 1 + P lj {AQ ■ Vu) 1 , 

which can be rewritten as 

PA 2 u = A(PP ■ Vu) + div {(PE - 2VPD) • Vu) 

+ (APP> + PG + PAtt - VP • E) ■ Vu. 
Substituting (13.41) into (13. 3ft . after some rearranging, yields 

A(PAw) = A(PP • Vu) + 2div 2 (VP <g> Vu) 
(3.5) + div ((PE - 2VPD - 2V 2 P • Vu) • Vu - APVu) 

+ (APP + PG + PAO + V(AP) - VP • E) ■ Vu. 



Define 



(3-6) E% = P^E% - 2d a P l W^ - 2d a d p P a - 5 aP AP u 

H = APD + PG + PAVL + V(AP) - VP ■ E - *dA£P, 

where, as before, the Roman lowercase indices run from 1 to n, and the Greek indices from 
1 to 4; 5 a p denotes Kronecker's delta. 
With this notation, we rewrite (13.51) as 

A(PAw) = A(P • Vu) + 2div 2 (VP <8> Vu) 

(3.7) 

+ div(P • Vu) + H ■ Vu + *dA£ ■ PVu. 

We shall need the precise form and integrability properties of terms that appear in H. By 
Theorem 12. 1} 

*dA£P = A ((VP + P^)P _1 ) P, 

thus 

H = APD + PG + PAQ + V(AP) - VP ■ E - A ((VP + Pfi)P~ 1 ) P, 
and after some simple reductions we get 

H ii = A pu D i3 + pii G ij _ d a p il E% - ap'oJ' 



- 2(d a P u d a ^ 3 + d a dpP u d a (p- l ) ls P s i 

+ d a p il n lt p d a (p- l ) ts p s i + p il d a np p (p- 1 ) ts p sj 

- d^p il A(p- 1 ) ls p sj - p u n l ^A(p- 1 ) ts p sj , 
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or, in the simplified notation, without the jungle of indices, 



(3.8) 



H = APD + PG -VP • E - APtt 

- 2(VP • Vfi + (V 2 P • V(P~ 1 ))P 

+ (vpfi) ■ v(p _1 )p + P(vn ■ v(p~ 1 ))p 

- VPA(P" 1 )P - P^A(P- 1 )P 



Please bear in mind that in fact we shall use equations ( I3.6|) — (j3.8|> for each k, adding the 
subscript k to all letters H, D, E, G, P, ft, £, D, E. 

From now on we shall return to using the index k, where appropriate. 

The bounds for Dk, Ek, Gk and Qk together with the estimates for Pk and given by 
Theorem 12.11 imply that Hk G L 4 ^ 3 and 

(3.9) \\H k \\ L 4/3 <C = C( sup max (||£> fe || w i,2, ||P fe || L 2, ||G fc || L 4/ 3 , ||Q fe || w i, 2 ) J . 



To check this, one just uses Holder's inequality, and - when it is appropriate - the Sobolev 
imbedding W 1 ' 2 C L A in dimension 4; we leave the details to the reader. 

Remark. If the coefficients Dk = D(uk, Vrtfc), Ek = E{uk,Vuk,V 2 Uk), Gk = Giuk.Vuk) 
and Qk — &{uk, Vu/c) are given by the composition of fixed smooth functions with the Uk 
and their derivatives, and satisfy the growth conditions (jl.lljl . then a computation yields 



We do not rely on that particular estimate, though. 

Let H be defined analogously to Hk-, i.e. by formula f)3.8p j^ The convergence of all the 
terms in the right hand side of the formula for Hk is such that, up to a subsequence, we 
can assume 



To see this, we just use the elementary observation: if 1 + ^ = p,q,r > 1, and fk — >■ / in 
LP and gk — ^ g in L q , then ]k9k fg in L r , combining it with the imbedding W 2,2 C W 1 ' 4 " 
and Rellich-Kondrashov's compactness theorem. 

Moreover, by the estimates on Hk, we can once again choose a subsequence of Hk that 
is convergent weakly in L 4 / 3 . 



k 




<||n fc || 3 ^, 2 + ||V« fc ||l 4 . 



(3.10) 



H k ^H in L q for aU 1 < q < -. 



'This time the subscript k in (|3.8p is really omitted, not just for the sake of simplicity! 
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We write out the weak formulation of (jl.4p in Bj, using its transformed form (13. 7p . and 
separating the terms into 'easy' (left hand side) and 'hard' (right hand side): the identity 

J P k Au k AiP - j(b k - Wu k )Aij 
(3.11) + J [2div(VP fc ® Vu k ) + (E k ■ Vu k 

(H k ■ Vu k )^ + J *dA£ k ■ {P k Vu k )^ 

holds for each smooth map ip compactly supported in B = Bj. (Since £ k is only of class 
W 2 ' 2 , the last term has to be interpreted using one integration by parts.) 

3.3. Convergence of ( 13. lip . We consider the left hand side and the right hand side 
separately. (The key difficulty is to prove that right hand sides converge to the appropriate 
limit). 

3.3.1. Convergence of the left hand side of (13. lip . By assumption, u k — u in W 2 ' 2 , thus, 
after passing to a subsequence, we may assume that 

u k — > u in L s for all 1 < s < oo, 

Vwfc — 51 Vw in L 4 , 

Vu fc V« in L p for all 1 < p < 4, 

V 2 u k -± V 2 u in L 2 . 

Taking into account the bounds for D k , E k , G k and fl k we can also assume that (again, 
after passing to a subsequence) there exist D, E, G and Q such that 

D k — ^ D and tt k Q in W 1 ' 2 , 
E k ^ E in LP for all 1 < p < 2, G k -± G in L s for aU 1 < s < 4/3. 

Similarly, Theorem 12.11 gives uniform estimates on P k and which allow us to assume 
(after passing to subsequences) that 

(P fc ,6)^(P,0 in^ 2 ' 2 , 

and, by the Rellich-Kondrashov's compactness theorem, Vpt — > VP and — > V£ 
strongly in L s for all 1 < s < 4, and P k — > P in L s for all 1 < s < oo. The limits P and £ 
of P k and £fc satisfy the claim of Corollary 12.31 on B = Bj. In particular, 

dPP- 1 + POP" 1 = *d£. 

By Holder's inequality, it follows from all these convergence assumptions that VP k ®Vu k — > 
VP ® V« in P', 

. Vwfc = PfcP fc ■ Vw fc PP • Vu = P • Vtt in P', 

and finally, 

P fe • Vw fe -> P • Vw in P' 
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(to check this, see the (13. 6ft for relation between E and D, E, and P). 

Thus, up to passing to a subsequence, for each fixed test map ip the left hand side of 
(13. lip converges to 



J PAuAtp - J(D ■ Vw)A^ + J 2div(VP ® Vtt) + E ■ Vu 



3.3.2. Convergence of the right hand side of (13. lip . Now let us concentrate on the conver- 
gence of the right hand side of ( 13. lip , i.e. of 

(3.12) J (H k ■ Vu k )*p + J *dA& ■ (P k Vu k )ij. 

This is the heart of the matter. A priori, by Holder's inequality, H k ■ Vu k is bounded just 
in L 1 , and there is no simple means of passing to the limit, as neither H k — )■ H in L 4//3 , 
nor Vufc — > Vw in L 4 ; the convergence in both cases is weak. The second term presents 
a similar problem. To circumvent this difficulty, we shall study the convergence of (I3.12p 
using the concentration-compactness method of P.-L. Lions. To deal with the first term in 
(I3.12p . we exploit the fact that u k has second order derivatives and use Sobolev inequality. 
To cope with the second term, we proceed in a similar way, employing also the equations 
satisfied by P and £, and integration by parts. 
We define the auxiliary distributions T k by 

(3.13) (T k ,tff) = J(H k - Vu k )*p + J *dA£ k ■ (P k Vu k )^, ^ e C™(B,R n ). 
Coordinatewise, we write T k = (T k , . . . , T£), with T % k G V(R n ) given by 

(Tl <p)= j {Hi ■ Vu k )<p + J *dAC k ■ (P k Vu k )<p, <\> e C™{B). 

We define the distribution T by a formula analogous to ( I3.13p . omitting the index k every- 
where. 

Let us recall again that, for % = 1, . . . , n, 

W = (2#), 
*dC = 

with j — 1, . . . , n and a = 1, . . . , 4. In view of that we have 

(pH^Vu = (f)H l Jd a u j 
<P * dAC ■ PVu = <f)(* dAi)l P jk d a u l . 

Lemma 3.1. There exists a subsequence T k > out of T k that 

T k , ^T l + J2 a u 5 *u in V'(R n ), i = 1, . . . , n, 
leJ 

where J is (at most) countable, an G K, Xn G B 4 and Y^iej l a z«| 4 ^ 5 < 00 f or i = 1, ■ ■ ■ ,n. 
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Proof. Our aim is to prove that there exist non-negative, uniformly bounded Borel measures 
/ifc such that for any <fi e C™(B) and i — 1, . . . , n we have 

(3.14) \(Ti-T\tf)\<(J (frU^' +o(l) as A; ^ oo. 

This will allow us to use the method of concentrated compactness. 
We have 

(7* - T\ 5 ) = y 5 (iZj • Vu k - H* • Vu) 

( 3 - 15 ) + y 5 (*dA£* • P fc Vw fc - *dAf • PVu) 

= A + B. 

Each of the integrals shall be dealt with separately. 
Estimate of A: we split this integral into two, 

A = J cj) b H l k ■ {Vu k - Vu) + J (t> 5 {Hi - H' 1 ) -Vu=I + II. 

For I we have 

/ = J <P 3 HI ■ V (0 2 K - «)) - 2 y 4 ^-(V0®K-n)). 

We have assumed that — >• u in L s for 1 < s < oo and that \7u k — > Vu in L p for 
1 < p < 4, moreover, we know that the H k are uniformly bounded in L 4//3 . By Holder's 
and Sobolev's inequalities 

|/| < y \<f)\ 3 \H k \\V((P 2 (u k -u))\+o(l) 

< (y^i^i 473 ) 374 (y iv(0 2 K-n))i 4 y /4 +o(i) 



(3.16) < (J 4 |^.| 4/3 ) 3/4 (y |V 2 (0 2 K-n))| 2 y /2 + o(l) 

\ 3/4 / - \ 1/2 

4 |iJ fc | 4 / 3 J (/ 4 |V 2 K-n)| 2 J +o(l) 



< 

no 



< 



|^| 4/3 + |V 2 n fc | 2 + |V 2 n| 2 )) 



5/4 

+ 0(1). 



With 77 we proceed in a similar way: 

11 = J ^{Hl - H*) ■ V(0 2 u) - 2 y 4 (/7* - W) ■ (V0 ® u) 



(3.17) 



II a + // 6 
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We estimate II a in the same way as we did with /: 
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Wa\< (|0 4 (l^| 4/3 + l^| 4/3 +|V 2 H 2 )) 



5/4 



o(V 



The integral IL converges (up to choosing a subsequence) to zero, since u G L s for any 
s > 1 and H k ^> H in L q for q < |. 

Estimate of 5: this integral is the sum of three terms, 



5 



(3.18) 



y 5 (*dA£ ■ P k Wu k - *dAC • PVu) 

* c/A^ • P k V{u k -u)+ [ 5 * dA£* • (P fc - P)Vu 



+ / 5 (*dA^j - *dAf ) ■ PV 



/// + IV + K. 



Since, thanks to the boundary conditions on A£ 



J *dAC k ■ V {<fP k (u k - «)) = 0, 



we obtain 



(3.19) 



</> 4 *dA£*- (V0®P fc K-n)) 
i 5 *dA^-(VP0K-u) 

P/a + ///ft. 



Integrating by parts, taking into account the convergence of u k to u in any L s and arbitrary 
integrability of P k we obtain 



\IIL 



(3.20) 



< 



Aei Ad^p^uk-u)) Ad 

J <p 4 \V<P\\V%\\VP k \\u k -u\ + J 4 |V0||V 2 6||P fc ||V Mfc -V 



o(V 
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For Illf, we get, integrating by parts, 



\IIIh 



AC k dP k Ad(<j) 5 ( Uk - U )) 



(3.21) 



< 



< 



< 



*A£dP k A d(u k - u) 



l \V4>\\V%\\VP k \\u k -u\ 



{<p 2 AC k ) (<f)dP k )A(<p 2 d(u k -u)) 
(0 2 A^) {(p dP k ) A d (4> 2 (u k — u)) 



oil) 



o(l). 



As before, using Holder's and Sobolev's inequalities, we estimate the last integral by the 
product 



1/2 / f \ 1/4 

l4|V72 f \1\ (l J.4IV7 D |4 



4 lvp fe l 



V z {(j>\u k -u))\ 



1/2 



This yields 



(3-22) II h < (| 0^ (|V 2 a| 2 + |VP fc | 4 + |VV| 2 + |V 2 w| 2 )^ 



5/4 



0(1). 



The remaining integrals in B - the terms IV and V - are estimated in much the same 
way. Integrating by parts, and then dealing as in the proof of f)3.22p . we obtain 



(3.23) 



\IV\ 



5 * dAC k ■ (P k - P)Vu 
Aad(<P 5 (Pk-P)))Adu 



< 



^ 2 AC k d(<P(P k - P)) A d{ft 



+ AJ ^\V<P\\V^ k \\P k -P\\Vu\ + o(l) 
< (/ 4 (I V 2 6| 2 + |VP fc | 4 + |VP| 4 + |V 2 w| 2 )^ 



5/4 



0(1). 
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Similarly, 



V\ 



(3.24) 



< 



A(a-C)d(<f> 5 P)Adu 

<j> 5 A(a-e)dPAdu 



< (7V (|V 2 6| 2 + |V 2 £| 2 + | VP| 4 + |V 2 w| 2 )^ 



A(C k - C ) A d(j) A Pdu 

5/4 



0(1). 



(in the calculations above we use Holder's and Sobolev's inequalities together with the fact 
that V 2 £ fc V 2 £ in L 2 ). 

Altogether we obtain the estimate 

\B\< (J 4 / fc ) ' +o(l), 

where f k := |V 2 ^-| 2 + |V 2 £| 2 + | VP fe | 4 + |VP| 4 + | V 2 w fc | 2 + |V 2 w| 2 E L 1 . Putting together 
the estimates for A and B we get 

5/4 

(3.25) 



{n-T\<f)\<(^J ^d^ +o(i) 



with 

(3.26) ti k = \H k \ A / 3 + \H\ A / 3 + f k . 

Note that by our assumptions the fi k are uniformly bounded in L 1 . Passing to the limit in 
the space of measures we obtain 



T k -T -> du, 



and 



'du 



< 



d[i k — > du, 

\ 5/4 

b*dn J . 



Now the claim of Lemma 13.11 follows directly from the concentration-compactness lemma 
of P.-L. Lions, cf. [SJ Lemma 1.2, p. 161]. □ 



Passing to the limit in (13. lip , we obtain 
PAuAip - / (D ■ Vu)Aip 



(3.27) 



+ 



2div(VP <8> Vu) + (E ■ Vu) 



J (H ■ Vu)i) + J *dA£ ■ (PVm)^ + (Sj, i)) 
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for each smooth test map ip compactly supported in Bj, with H, D and E given by ( 13. 6p . 
and (P, £) related to Q = UmQk via Uhlenbeck's Theorem 12.11 The singular distribution 
Sj G T>' is a countable series of Dirac delta measures with vector- valued coefficients, Sj = 
^2i a l$xf If follows from Lemma [3. II that the series ^2 \ae\i converges. 

3.4. Removing the singularities. To show that the limit u of u k satisfies the limiting 

system not just in Bj \ A, where A = {xe: I G J} is countable, but in fact in the whole 
Bj, we rely on the fact that W 2 ' 2 (M. i ) contains unbounded functions. More precisely, the 
following holds. 

Lemma 3.2. There exists a sequence of functions G C^°(1R 4 ) such that < < 1, 
$fc = 1 on B(0, Tfc) for some r k > 0, 

(3.28) $ fc = on W n \ B(0, R k ), R k -t as k -> oo, 
and 

(3.29) ||$ fc ||vi/2.2 -> as k ->■ oo. 



For the reader's convenience, a proof of this lemma is sketched in the appendix. 

Fix £o G J and assume for the sake of simplicity that X£ = G Bj. Testing each 
equation % — 1, 2, . . . , n of system (I3.27P with vjjk = (±$fc, • • • , ±$fc), the signs being equal 
to the signs of coordinates of the coefficient ae in Sj, and keeping in mind that 

PAu, D ■ Vu, VP ® Vm, A£ g L 2 , H-VueL 1 , d((PVu)4) k ) e L 2 , 
we easily obtain 

o(l) = o(l) + (Sj,ijj k ) = o(l) + | a £o 1 1 + ^{ae^kjxi)) , k ->■ oo. 

By (I3.28P and Lebesgue's dominated convergence theorem, the sum X^£ a bove tends to 
as /c — 7- oo. Thus, upon passing to the limit, we obtain = for i = £ . It follows that 
Sj = 0. 

Thus, relying on the above and selecting, via the standard diagonal procedure, a sub- 
sequence k' — > oo such that u' k — > u (and all the coefficient functions converge in their 
appropriate spaces) in all good balls Bj simultaneously, we check that 

(3.30) A 2 u = A(D ■ Vm) + div(E ■ Vm) + G ■ Vu + AQ ■ Vu in B 4 \ A u 

where A\ = {x\, . . . ,x s }, s < N, denotes the finite set of bad points, cf. (13. 2p . To see that 
in fact (I3.30P holds in the whole ball B A , we pick an arbitrary test function ip G C^°(M A ) 
and write it as the sum 

<P(?) = P( x )®k(x - Xj) + (p(x) f 1 - ^ $ fc (x - Xj) J =: ip ,k{x) + <pi t k(z). 

j=l V j=l ' 



WEAK COMPACTNESS FOR 4-TH ORDER SYSTEMS 17 

Then, tpi^ — > <p in W 2,2 as k — > oo. Moreover, supp^x^ <e B 4 \ A\. It follows easily that 
the weak form of ( I3.30P holds not just for test functions tpx^ G C^°(B 4 \ Ai), but also for 
an arbitrary ip G C^°(B 4 ). 

The whole proof of Theorem 11.11 is complete now. 

4. An EXAMPLE 

The following example shows that a system of type ( II. (jp satisfying assumptions (jl.3p 
may have solutions that are in C 1 \ C 2 . 

Set 0(r) = loglog(e — logr) for r > 0. 

Lemma 4.1. TTje function f(x) = <j>(\x\) = log(log(e — log(|x|))) is in W 2 ' 2 (M 4 ). 

Proof. The function is obviously in L 2 (B 4 ), therefore (e.g. by Gagliardo-Nirenberg's in- 
equality) it is enough to check that V 2 / G L 2 (B 4 ). This is done by an elementary compu- 
tation which, for the sake of completeness, is sketched in the appendix. □ 

Lemma 4.2. The functions t<p'(t), t 2 (j)"(t) and t 3 4>"'(t) are bounded on (0, 1]. 

Proof. A computation shows that 



\t 2 <p"{t)\ 

|tV"(t)l 



(e — logt) log(e — logt) ' 
— 1 + (e — 1 — logi) log(e — logt) 

(e — log t) 2 log 2 (e — log t) 
-2 + 3 (-1 + e - log(t)) log(e - log(t)) 
(e-log(t)) 3 log(e-log(t)) 3 
(2 - 3 e + 2 e 2 + (3 - 4 e) log(t) + 2 log(t) 2 ) log(e - log(t)) 5 



(e-log(t)) 3 log(e-log(t)) 3 
Clearly, the right hand sides converge to as t — > 0. The lemma follows. □ 

Example 4.3. Define w : B 4 -> R, 

(4.1) w(x) = \x\ 2 sin 0(|x|). 

One easily checks, using Lemma [4. 2 1 that Vu> and D 2 w are bounded on B 4 . However, 

d 2 w 

-(0,0, 0,t) = 2sin0(t) +3t<f>'(t) cos cj)(t) 



dx\ 

does not have a limit for t — > 0, and thus w is not a function of class C 2 . Finally, 

|V 3 ^(x)| < , . G L 4 (B 4 ). 

|x|(e — log \x\) log(e — log \x\) 

Set C > to be the bound on |Viu| on B and consider 

4 

(4.2) v(x) = w(x) + 2C^Xi , a(x) = Av\Vv\~ 2 . 

i=i 
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By definition, |Vi> | > C > on B 4 , and \Vv\ and \D 2 v\ are bounded on B 4 . Thus, a(x) is 
bounded on B 4 . 

We also have \D 3 v\ = \D 3 w\ G L 4 (B 4 ), and 



iVal < 



\VAv\ 




2 + l 


Av\ 


\Vv\ 


\D 2 v\ 




Vv 


4 



< 1 + \D 3 v\ e L 4 



Thus, a(x) G W 1A (M A ) n L°°(B 4 ). Consider the equation 

(4.3) A(Au(x)) = A (a(x)\Vu\ 2 ) 
Clearly this is an equation of the type 

(4.4) A 2 u = A(D (u, Vw) • Vw) + other terms (with zero coefficients), 

where D = a(i)V?i satisfies 

• D E W 1,2 whenever u G W 2 ' 2 , 

• |-D( 3: )| < |l a IU°°|Vu(x)|, 

• ||V-D||l2 < || Va|| L 4||VM||L 4 + ||a|U°°||-D 2 M|U 2 - 

It is also obvious, by the very definition of a, that v solves (14. 3p . Therefore (I4.4p . under 
the conditions on D listed above, admits non-smooth solutions. 

Note, however, that the function D given above does not satisfy the pointwise estimate 
for VD given in (II. lip . 



Auxiliary lemmata 

An unbounded function in W 2 ' 2 (M A ). As explained in Section 4, to show that f(x) = 
log (log(e — log |a;|)) is in W 2,2 on the unit ball B 4 , it is enough to check that D 2 f G L 2 . 
We have (Mathematica) for x G B 4 : 

dl x J(x) 



x i x i (l + (l-2e + 21og 


x\) log(e - 


- log 


A)) 




X 


4 (e - log 


x\) 2 lo{ 


? 2 (e-log 


x\) 



< 



•2/ ^ OC j 



|x| 4 (e — log \x\) log(e — log \x\) 

< ! . 

~ |x| 2 (— log |x|) log (— log |x|) 



Integrating the above squared over B 4 , in polar coordinates, amounts to calculating 



dr 



dt 



io r log 2 r log 2 (— logr) j\ t log 2 1 log 2 log t ' 
with the latter integral convergent, since the series 

1 



E 



n log 2 n log 2 log n 



is convergent, by Cauchy's condensation test. 
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Proof of Lemma 13.21 Let / £ W 2 ' 2 be the function introduced in Lemma 14. 11 Fix a 
function ( £ C°°(R) with < ( < 1, C = on [1, +oo) and C = 1 on (-oo, |). Set 

(A.l) $ fc (x) = C(A;+l-/(x)), x£B 4 , fe = 1,2,... 

Clearly, $ £ W 2,2 (M 4 ). Moreover, $ k (x) = if f(x) < k. The equality $ k (x) = 1 holds 
whenever f(x) > k + |, i.e. on a neighborhood of 0. Since / > on B 4 and / is smooth 
except at 0, $fc is smooth on B 4 , and can be extended by to the whole space M 4 . Finally, 

|V**(*)| < |V/(x)| -X {f>kV 
|V 2 $ fe (x)|< (|V 2 /(*)| + |V/(z)| 2 ) - X{f>k} . 

Since / £ W 2,2 fl V^ 1,4 is radial and / — > +oo at 0, we obtain ||<&fc||w2,2 — > as k — > oo. 
For x y£ 0, we simply have &k{x) = for all k sufficiently large. This completes the proof. 
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